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Abstract 

The sector of zero Zjv-charge is studied for the ferromagnetic (FM) and an- 
tiferromagnetic (AFM) version of the Zn x Z2 invariant Fateev-Zamolodchikov 
quantum spin chain. We conjecture that the relevant Bethe ansatz equations 
should admit, beside the usual string-like solutions, exceptional multiplets, 
and a number of non-physical solutions. Once the physical ones are identified, 
we show how to get completeness and the gapless excitation spectrum. The 
central charge is computed from the specific heat and found to be c = 
(FM) and c = 1 (AFM). 
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1 Introduction 



It is of interest to study the Fateev-Zamolodchikov spin chain 

M N-1 -1 

where the operators Xk and Zk act on the A^*^-dimensional vector space spanned 
by the basis n= \ni,n2, ■ ■ ■ um >, Uk = 0,1 . . . N — 1 

Xk\ni . . .Uk . . .riM > = \ni . . .Uk + I . . .Um > modA^ 

Zk\ni . . .Uk . ■ -Um > = u!^''\ni . . .Uk . . .Um > w = exp{2TTi/N) 

in its ferromagnetic (i.e. H itself) and its antiferromagnetic version (i.e. —H). 
The Hamiltonian (0) is exactly integrable since it can be derived from the family of 
commuting transfer matrices T{u) of an integrable 2-dimensional spin model 

[T{u),T{u')] = yu,u'eC 

N-1 1 

^1 sm(n7r/iV) 

As it was conjectured in the paper 0, and subsequently proved through the calcu- 
lations of the critical exponents [Q, the model is critical and, in the scaling limit, it 
gives a Zj^ invariant conformal field theory with parafermion currents. 

While this leaves little to be discovered about the ferromagnetic (FM) version of 
(|^), little is still known about the antiferromagnetic version (AFM). The case A^ = 3, 
corresponding to the three-state critical Potts chain, has been studied in great detail, 
leading to the calculation of the excitation spectrum 0, the central charge and the 
spectrum of conformal dimensions of the corresponding field theory ^ and finally 
the characters of the relevant representations of the Virasoro algebra as well as the 
full modular invariant partition function 0. On one hand, these results have led 
to recognize that the A^ = 3 AFM spin chain is critical and in the scaling limit it 
describes, rather surprisingly, a conformal field theory with Z^ parafermions. On 
the other hand, even the results for A^ = 3 FM have proven fruitful because they 



have provided a new way, directly related to lattice models, to express characters of 
the representations of the Virasoro algebra and affine Lie algebras 0. 

In this paper we extend some of these results to arbitrary (but odd) N. We 
summarize here what is known about the exact diagonalization of (|I]) 0. By means 
of the analitic Bethe ansatz, the whole spectrum of (|l]) can be expressed in terms 
of a set of variables {A^}, related in a simple way to the zeroes of the transfer 
matrix eigenvalues of the two dimensional model. The set {A^} obeys a system of 
trascendental equations whose appearance is commonplace in integrable models 

sinh(Aj - Afc + i-f) ^ ^.^^^ sinh(Aj + ^^7) , s 

^}^sm\i{Xj- Xk-i^) ^ ' ^sinh(Aj -zs7)^ 

(iV-l)7r 1 
7 = — — - — s 



2N 2{N-l) 
L = {N~l)M-2Q g = 0, l...(A^- l)/2 

Here Q is the Zj^ charge, defined modA^, corresponding to the eigenvalue of the 
conserved operator 

M 

exp(227rg/iV) = X{Xk 

k=l 

Sectors of charge Q and N — Q are mapped into each other by the (conserved) charge 
conjugation operator 

C\ni, n2, . . . um >= \N - Til, N - n2, . . . N - um > 

so that the symmetry group of (|l|) is Z/^r x Z2. Energy and momentum of each state 
are related to a solution of (0) 

L {N-l)/2 

E = ^cot(zAfc + — ) -2M J2 cot{nk/N) (3) 
k=i k=i 

^ sinh(Afc + ^) 

exp(^P) = 11— TTT 1^ (4) 

fc=i smlil^Afc - j^) 

In the following sections we consider the sector Q = and show that the excitation 
spectrum is massless and made up of one species of "quasiparticle" (FM) and 



3 



species of "quasiparticles" (AFM) with linear dispersion relations at small momenta 

FM e{p) ^ Np (5) 
AFM e{p) ^ j^p (6) 

We also determine a set of completeness rules that allow to classify all physical 
solutions of (0), i.e those solutions that actually correspond to an eigenstate of (Q), 
and show how to count them to obtain the correct dimension of the sector Q = 
in the vector space spanned by the basis n. Finally, the set of nonlinear integral 
equations |jlO[ that describe the thermodynamics of (|l|), are used to determine the 



low temperature asymptotics of the specific heat, related to the central charge of 



the conformal field theory by |11 



3v 

where v is the group velocity of massless excitation, from which we find 

A^ - 1 

FM c = 2- (7) 

AFM c = 1 (8) 

The value (|^ is of course the one predicted in 0. 

2 Strings and multiplets 

The traditional hypothesis on the solutions of (0) is that, in the limit M —>■ oo, they 
group into strings 



A!:/) = A(-) + |(n + l-20 + ^^^^^ 1 = 1,2. ..n (9) 
where A^"'*"^ is the real center of the string, n its length, and v = ±1 its parity. It's 



been long known |T2|, and the numerical analysis of the case A^ = 3 confirms it |]T3 
that complex pairs with imaginary part different from the one given in can also 
appear. We formulate the following conjecture: for fixed A^ = 2p + 1, the solutions 
(roots) of (0) fall into three classes 



1. 1-strings with both parities: f = ±1 

2. Even length strings: {n,v), n = 2,4, . . . N — 1, v = ±1 

3. p multiplets, to be denoted (M, m), of length 4m + 2, m = 0, 1, . . .p — 1 

A!,^-) = Af-)±.(^ + i^) l = -m,-m + l,...m-l,m (10) 

with A^*^'™^ real. So, for = 3 = 1) we have only a pair (m = 0), for = 5 
(p = 2) we have a pair (m = 0) and a sextet (m = 1), etc.. This conjecture is 
partly motivated by a numerical diagonalization of the transfer matrix and partly 
warranted a posteriori by the fact that it produces the correct counting of states 
(see Section 4). To diagonalize T{u) numerically, one fixes the spectral parameter u 
at some conveniently chosen complex value, finds the eigenvectors numerically and 
then applies T{u) to them. The eigenvalues are then polynomials in e^*" and it is 



easy to locate numerically their zeroes ||T3[| . 



Equations for the real centers are obtained by taking the product of (0) over 
members of a string (multiplet), so that all factors are in the form 

_ sinh(^a + A + z(l-t;)f) 
^ '"'""^ sinh(2a-A + z(l-^;)f) 
XeR ae (-7r/2,7r/2] 

In (^) strings have been assigned their limiting "perfect" value, and one may have 
a = or a = tt/2 with G{X,0,v) = —v and G{X,7!'/2,v) = v. We group these 
exceptional factors to yield an overall ± sign, and then take the Log choosing the 
branch 

iLogG'(A, a, v) = 0(A, a, v) = 2v arctan(cot(a)'' tanh(A)) (11) 
a G (-vr/2,0) U (0,7r/2) 

so that (j){X,a,v) is a monotonical continuous function with values in (— 7r,7r). The 
new equations are written with the help of Z-functions 

^.(Ai^-)) = J^Vm (12) 



1 

2^ 



1_, 

I ■ 

k p=l 



27rM ^ ^ ^'^ ^ ^ 



(13) 



Here j, k denote the type of string (multiplet), is their number, l'^^ are integral 
or half-odd, and the functions in (|TB|) are defined 



TT 



t(M)(A) = 2<P{X,—,v) 



t{M,m)W 



J22<f){X,'^{n + 2s + l-2l),v) 
1=1 

, , ,^71 2m + 1 
^ 2e0(A, - + e- 

e=±l ^ 



4A^ ' 



(14) 
(15) 

(16) 



0(n,t.),(n',«,)(A) = 0(A, -(ra + n'),t;u;) + 0(A, -|n-n'|,ww) 



min(?i,n') — 1 



+ 2<P{\,]-{\n-n' + 2l),vw) 

1=1 ^ 



(including the case n or n' =1) 

e 



(i,«),(M,m)(A) = y [e0(A, 7 

e=±l 



4 2iV'") + ^'^(^'4 



TT e7r(m + 1) 



2N 



(17) 



^)] (18) 



2m 



0(n,«),(M,m)(A) 



e=±l ^ 



©(«,■«), (M,m) (A) 
e=±l ^ 



E E20(A,^vr--^(n-2m-l + 2/),t;) + 



_ 2m - 1), t;) + 0(A, - ^(n + 2m + 1), v)\ 

E E20(A,i^7r--^(2m + 2/ + l-n),t;) + 
_^(2m - n + 1), t;) + 0(A, - ^(2m + n + 1), v)\ 



(19) 
(n > 2m) 

(20) 
in < 2m) 



,m),{M,m' )(A) 



2 E [<^(A, - ^l'^ - ^'1, +) + 0(A, - ^(m + m' + 1), +)] 



2min(m,m') ]^ _|_ g ^ 

+ 4E E '/'(A, ^vr--^(|m-m'| +/),+) 

e=±l 1=1 ^ 



(21) 



The notation in (|l^pi| ) has been used for the sake of brevity but it is understood, in 



agreememt with the remarks before ([TT|) , that in (/)(A, a, v) a is shifted by periodicity 
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to be in (— f,f] and 0(A, 0,f) = 0(A,7r/2,f) = 0. These exceptional values of a 
only change the oddness of l'^^^ in (|12D and keeping track of them one finds 

/(!'+) = integer(half - odd) if 1 + M + M(i,+) + Mjv-i,(_)p) = even(odd) 

= integer (half - odd) if 1 + M + M(i _) + M(^_i^(_)p+i) = even(odd) 

j(7v-i,(-)f+i) _ integer(half - odd) if 1 + M(i,_) + M(^_i,(_)p+i) = even(odd) 

j{N-U-r) ^ integer(half - odd) if 1 + M(i,+) + M^n-U-)") = even(odd) (22) 

= integer(half - odd) if 1 + M(„,+) = even(odd) 

j(n -) ^ integer(half - odd) if 1 + M(„ _) = even(odd) 

j(M,m) _ g^j^g^yg integer 

Energy (^ and momentum (^) become functions of the string centers. Real part of 
the energy and momentum are simply related to the t-functions in (|1^|16]) 

Ree,(A) = h'{X) Vj 

p,{X) = -U, VjV(1,+) (mod27r) 
P(i,+)(A) = -it(i,+)(A) + 7r (mod27r) 

but even after summation over string (multiplet) members, the energies retain an 
imaginary part 



sinh(2A) 

sh(2A) - cos(^ 

including (1, ±)-strings and 



Ime(„,,) g cosh(2A) - cos(^7r(n + l-2k) + ^vr - ^) ^^^^ 



_ ™ sinh(2A) sinh(2A) 

" " ,£r.^cosh(2A) + sin(^ + cosh(2A) + sin(^? ^''^ 

It will be shown, at least in the thermodynamic limit, how these imaginary parts 
can be removed assuming suitable correlations between the quantum numbers l''^^ 
in 
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3 Ground state and excitations. Correlation of 
rapidities. 



The survival of an imaginary part in the bare energies ( P3| - |2^ signals that correla- 
tions between the rapidities {A^^^} must exsist in order to ensure reality of the total 
energy. In fact, the detailed investigation carried out in the simplest = 3 case 
T3| clearly shows that (|I2|) contain several spurious solutions and only a subset of 
the possible choices of quantum numbers {I^^ in (|T^ reproduces the correct phys- 
ical solutions, i.e. the spectrum of H. This peculiarity, together with the massive 
appearance of non-string multiplets, is probably a consequence of the "unorthodox" 
approach used here to diagonalize the transfer matrix |p. In fact, the unknowns of 
(0) are, up to a change of variables, the zeroes of the transfer matrix eigenvalues 
themselves, rather than the zeroes of an auxiliary Q-matrix satisfying a "T — Q 
recursion relation" |T^, |T3[. Therefore some details of the conventional method of 



dealing with (|T2D must be modified. In particular, densities of rapidities and holes 



that describe the solutions of ([T2|) in the thermodynamic limit are expected to be 
related by a set of constraints ^ and this should be taken into account if one were 
searching for a minimum of the free energy functional of the gas of rapidities and 
holes . 



Instead of assuming from the start a set of selection rules on {ij^^, and conse- 
quently on the rapidities {A^-*}, we prefer to make a working assumption, essentially 
based on small chain observations, on the ground state in the infinite chain limit. 
We will check the validity of the assumption by performing small variations around 
the ground state configuration to ensure that it is indeed a minimum, at least lo- 
cally, of the energy. The case = 3 will provide an heuristic guidance. Since 
(A^ — 1, ±) and (1, ±) strings play a preminent role in the following, we adopt the 
shortened, A^-dependent notation: (1, +) = (a), (1, -) = (6), (A^ - 1, = (c), 

{N-l,{-y) = {d). 

The FM ground state is a band of (A^ — l,^—}^^^) = (c) strings M, whose 



jCAf i,( fQj^j-Q closely packed sequence symmetric around 0, and whose centers 
fill the real axis with density 

solution of ( here * means convolution) [|T^ 

P°(A) = -^:(A) = + ^Q'cc * P^c\^) (26) 

The (real) density of ground state energy is 

"° = i^oo aI = / ^^/'°(^)^^(^) - 2 E cot(7rfc/iV) 

fc=i 

In ( PB| ) we have taken into account that Zc is a decreasing function of A. Observables 
in the excited states are computed from dressing equations 0, ^ that incorporate 



the effect of the backflow in the ground state distribution (psf ) when holes are pinched 
in it and other types of roots are added. For a state described by a density of c- 
strings pc containing a finite set of M^'^'' holes {A^'^'^''}, and by a finite set {A^'^''} of 
other roots, we have 

1 Mi'") 

Z'^iX) = -a,(A) = -pe(A) - ^ E ^(A - ^f'""^) 



From the definition of Z-function 



cTcix) - ^e;,*o-,(A) = -^aA) 



1 m''^' -I Mk 



i: e;.(A - A--) + — E E eU(A - A-) (27) 

where the sum X^fc^-^c is taken over all roots other than c. Even without solving ( pTj) 
explicitly, one sees that the energy 

= EEefc(A?^)= / rfAa,(A)e,(A) 

fc /3=1 

- E^c(A(^"^) + EEe^(A?^) 

/3=1 fc^c/3=l 



can be written as 



Mi'"' Mfe 



E = Eo-Y. ^ciX^) + E E V) (28) 

/3=1 kjtcl3=l 

where Eq is the ground state energy and the dressed energies are defined as solutions 
of 

ek - ^Q'k,c *ec = ek (29) 

In (^) k runs over all strings and multiplets. The solution of (^) is easily found 
by Fourier transform method 

R«.(A) = -Ree4X) = -^^ (30) 

Imea(A) = Imec(A) = -A^tanh(2A^A) (31) 

all other being identically zero. The simplest way to cancel the imaginary part 
in (|28|) is to assume that pairwise 

{Xf'^^} = {Xf} p = l,2...M, (32) 



Then (^) is non-negative, thereby proving that ( P5| ) is indeed the ground state 
distribution (strictly speaking, it has been shown that it is a local minimum of the 
energy functional). It would seem that, due to the vanishing of all other dressed 
energies, one could add strings (or multiplets) other than (a) without increasing E, 
but it will be proved in the next section that this is not the case. 

The calculation of momentum is identical. The ground state momentum vanishes 
since p^°^(A) is an even function of A and in an excited state 

Mfe Mh 

P = M.vr + ^ ^ vr,(Af ) - ^ ^cixf^'^) (33) 

fe^c/3=l /3=1 

where the dressed momenta vr^ solve 

- ^0fe,c *^c= Pk (34) 
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It is convenient to work with a subtracted (odd in A) bare momentum p defined by 
Pk = Pk^k ^ (a) and Pa = Pa + tc, which accounts for the term Mail in (^). Then 
7r[, = — 2Reefc which is easily integrated 

7r„(A) = -7rc(A) = -2arctan(tanh(iVA)) 

and TTk{X) = for all other strings and multiplets. From (|28|), ( |30D and ( |32D it is 
clear that the excitation spectrum is made up of one kind of quasiparticle appearing 
as a "bound state" (not to be understood as a physical bound state of elementary 
particles): (1, +) paired to a hole in the (A^ — 1, (— )^^^) band 

cosh(2iVA;3) 

Ma 

P = Ek - 4arctan(tanh(A^A/3))] 

/3=1 

The dispersion curve for each quasiparticle is 

e(p) = 2Arsin(|) p G (0, 27r) 

which yields (^ at small momenta. 

The AFM case has a richer spectrum. We proceed as in the FM case, assuming 
a ground state, namely a filled band of (a) and a filled band of (6) corresponding to 
a closely packed sequence of J*^"-' and f^^^ in (|T^. Consequently, their densities are 
simply the derivatives of the relevant (increasing) Z-functions 

pf(A) = z;(A) = i-t;.(A)-^E0;v*pr^ 



2tx ' 27r 
where j takes values in {a, 6}. The solution is 



l=a,b 



27r J smh(l^) cosh( '' 

sinh(^^(^^ 



' 27ry sinh(|^)cosh(^^l ^ ^ 

The ground state energy density is 

eo = - E / dXpf\X)ej{X) + 2^ E' cot{k7i/N) (37) 



{N-l)/2 

^ E 

j=a,b k=l 
11 



We perturbe these distributions by pinching two sets of holes {A^ } and {A^ } 
and adding extra roots of type k ^ a,b 

M 



Zj(A) = a,(A) = p,(A) + ^ E '^(^ - ^?"^) 



/3=1 

where j G {a,b}. The following steps are a straightforward generalization of the 
previous case. Densities of vacancies are found from the system of two integral 
equations 

l=a,b 

^^^-"'(=0,6/3=1 ^^^-^ " fc^a.fe /3=1 

where j G {a, 6}. Owing to (|38D, energy and momentum in the thermodynamic limit 
reduce to 



E = i?o - E E e.(ArO+ E E^d>^r) (39) 

j=a,6 /3=1 fc7^a,fe /3=1 

p = Po + M.TT- E E^.(A?"^)+ E E^fclAf) (40) 

j=a,6 /3=1 k^a,bl3=l 

provided that the dressed observables e and vr are defined as solutions of 

1 

2^ 



efc + 7^ E e',,,. * e, = -efc (41) 

^fc + ^ E 01.,,*^^ =^;^ (42) 

where k runs over all types of roots. The solution reveals that all imaginary parts 
are odd functions of A and Ime(Af,fc) = 0, Ime(„,+) = — Ime(„,_), Imea = Imec, Imeft = 
Ime^, suggesting that, again, a pairing of rapidities must be in effect. One is led to 
assume that M(„,+) = M(„,_) = M„, = M^^), Md = Mf ^ and 

{AjT'^^l = {A?'-^} /5 = 1,2,...M„ (43) 
{Xf} = {Xf^-)} /5 = 1,2...M, (44) 
{Xf^} = {Xp'^} P=1,2...M, (45) 
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As to the real parts we have 

N 



RetaiX) = -Reee(A) 

AT - 1 S rn^h(m) + nosriHaik 



2 cosh(2A^A) 



2Ar cosh(f^) cos(^^^±i^^ 
,=oCosh(MA) + cos(i^) 



TV 

Reeb(A) = -Reerf(A) = Ree,(A) + 



4iV cosh(|^)cosr ^"-^^'+A)- 



cosh(2A^A) 

2(N-1) 



Ree. ±)(A) = — — V ( o L 2<n<N-3 

m ^ cosh(|^)cos(^^ 



^^""'"^^"iV-l ■±^.cosh(MA) + eos(|^ 



Ree^ 

It can be checked that Ree^ with j G {c,d, {n,v), (M, m)} is positive definite, con- 
firming that the conjecture about the ground state is correct. Clearly differ- 
ent massless excitations are present: p = ^^^^ multiplets, plus "bound states" 
{(n, +), (n, — )} correlated through (^) and finally the "bound states" {c, hole in a} 
and {d, hole in b} correlated through (^-J5). The equations for the momenta are 
easily integrated noticing that tt^ = 2Reefc 

sinh(f^) 

7ra(A) = -7rc(A) = - arctan(tanh(A^A)) - 2^ arctan( ^2j+i)nJ 

7r5(A) = -TTdiX) = 2arctan(tanh(A^A)) + 7r„(A) 

n/2 sinhC^) 
7r(.,^)(A) = 2^arctan( ^ ^^-Z'^ - ) 2<n<N-3 

j=l COS[ 2{N-1) ) 



7r{M,m)(A) = 2 arctan 



,sinh(-2^^ 
cosP^ 



j=-m ^^'^VaT-I/ 

These expressions for the dressed momenta are not sufficient to fix completely the 
dispersion curve for each quasiparticle since we have to deal with the contribution 
ttMq in (^Ol). This will be done in the next section after a rule on the string content 
of each physical solution of (ET) is derived. 
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4 Completeness. Correlation of quantum num- 
bers. 



The task of counting, let alone finding, all solutions of the system of equations 
(§) is very difficult. Nevertheless, in the framework of the string hypothesis, one 
can replace the original equations with ([T^) and hope that, to each set of distinct 
numbers {la^ there corresponds one and only one solution. Counting solutions 
becomes then a problem of combinatorics, and the result should equal the dimension 
of the vector space on which the Hamiltonian is defined. 

The sets {ij^^} cannot be chosen arbitrarily. An obvious constraint is that the 
total number of roots in the sector Q = must be 

N-3 p+l 

{N-1)M = M^+M,+iN-l)iM,+M,)+J2 E nM(„,,)+ ^ (4m+2)M(M,„^) (46) 

±1 rn=0 

even 

because L = {N — 1)M in (|^). Secondly, since the Z-functions are supposed to 
behave monotonically, the {I^^ must be chosen within the limits MZj(±oo). In 
principle, for a, say, increasing Z-function and a given string content of a state, 
one should determine the largest available vacancy Imax as the largest integer (or 
half-odd, according to (^)) strictly smaller than MZj{+oo) and count the integers 
(or half-odd) contained in [—Imax, Imax]- This gives the set of available vacancies for 
{la^, that we denote vac{j). It is a remarkable property of the sector Q = 0, and 
this is proven in Appendix B, that 

^ = ±(Z,(+oo)-Z,(-oo)) (47) 

where -l-(— ) is used if Zj is increasing (decreasing). 

An additional feature of (|12D is that, as it was extensively shown by the numerical 
study for the case = 3, only a subset of the possible choices {I^^, even when 
constrained by (^) and (^7]), gives an eigenvalue of H: the string content of an 
eigenvalue cannot be chosen arbitrarily and correlations are in effect between the 
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{ij^^- Motivated by these observations, and by ( ^5| - ^5D , we are led to assume that, 
for any state 

M(„,+) = M(„,_) = M„ (48) 
and denoting by {ij^'-'^} the empty vacancies (holes) and by Mj'^^ their number 

Mf ) = Me Mf ^ = Md (49) 

Since vac{j) = Mj + M^^\ imposing (|49| ) means 



Za{ + 00) - Za{-Oo) = — 

Zb{+oo) - Zb{-oo) = — 



Either one yields, when (^) and (|47|) are taken into account, the content rule 



Af-3 p-1 

(iV - 2)M, = NM, + 2{N - 1)M, + ^ 2nM„ + ^ (4m + 2)M(m,™) (50) 



n^2 ?ri=0 



which generalizes the content rule (3.1) of |Q, and it is in agreement with all nu- 
merical tests performed for N = 5, 7. It has to be understood as the statement that 
for any eigenvalue the numbers Mj must satisfy (^0]) as well as (^Gf). From (^) and 
(0) we choose to take Mb, M^, M„ and M(^M,m) as independent variables. We then 
find, after observing that all Z-functions except Za and Z^ are decreasing 

vac{a) = vac{c) = M + M^ + M^ 
vacib) = vac{d) = Mb + Ma 

vac{n, +) = vac{n, -) = ^^^(M^ + M^) + ^ ~ 2"^^^^"^' ^'))^n' (51) 
2,Ti^ n 2,771 -\- 1 

4:171 + 2 4777, + 2 71 

vac{M,m) = —-—{Mb + Md)+ E M^TTry - 2)M„ + E (4m + 2)(-— - - 1)M„ + 

^ n<2m ^ n>2m ^ 

^ , (4m + 2) (4m' + 2) , ,(4m + 2)2 ^ 

^ ^ 2(iV -2) ^ " ~ 2)M(M,„.') + i 2(N-2) ~ ~ ^)M(M,m) 
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The reduction by 1/2 of the vacancies of multiplets [|r^ has been adopted here too 
(see Appendix B). Of course, the number of vacancies must be, by definition, an 
integer. The content rule (|50D, implying that the right hand side must be divisible 
by — 2, guarantees that that is the case for vac{n) and fac(M, k). 

As to the correlations between {/^•'^}, from (|32|) and (|4^^5D , and noticing that 
fH6|) affects the oddness table fE2D, we are led to assume 



/(») = -/(M a = 1,2. 

m = AM) ^ = 19 



(52) 
(53) 



and, necessarily because of (pl|) 



jid) ^ _jih,b) c, = l,2...M, 



(54) 
(55) 



and from (^31) 



j{n,+) _ j{n-) 



(56) 



which generalize in a rather obvious way the selection rules of ||T^. It is shown in 
Appendix A that these conditions on the integers do imply the correlations ( p2D and 



The number of physical solutions of (|T2|) is therefore 



E 

Mt,Md,Mn,M(M,m) 



'vac{M^)\ f vac{Mb)\ ^'^ f vac{n)\ iJ (vac{M,m) 



Mr 



n 

even 



Mr,. 



n 

m=0 



M, 



{M,m) 



where 



M^ = M 



N-3 p-l 

-T7^{NMd + 2M6 + 5^ 2nM„ + ^ (4m + 2)M(m,™)) 

^ n=2 m=0 



(57) 



and the sum is subjected to the constraints 



1. the right hand side of (|i6| ) must be divisible by A^ — 2 



2. > 
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3. vac{j) > Mj \fj 

The result of (0) is expected to be N^''^~^, the dimension of the sector Q = 0. 
This has been checked numerically for =7, 9 and it will be proven analitically for 



= 5. The proof for = 3 has already appeared in |^ . 

It should be noticed that the number of vacancies is a piecewise continuous 
function of the imaginary parts of the roots belonging to a multiplet and it exhibits 
jumps at some special values. The conjecture (|TOD on the structure of the multiplets 
is based on the observation that it yields the right counting in (]^7\j. Therefore, the 
fact that (|57D give A^^^~^ should be regarded as a check of the various assumptions 
made so far. 

When = 5, (||) reduces to 

Ma = (5M, + 8Md + 4M2 + 2M(Af,o))/3 + 2M(Af,i) 
Hence, there must be an integer k, not necessarily positive, such that 

-Mb -Md + M2 + 2M(M,o) = -3/c 



which is used to eliminate M2 in favor of k. Then (|57|) becomes 



^ M + Mb + Md \/Mb + Md\ 

^[m -3Md-2Mb + 2M^M,o)+4:k-2M^M,i)) \ Md ) ^ ' 

/ Mb + Md + k \ / M(M,o) +2k\(2Mb + 2Md + M^m,i) \ 
\Mb + Md-2M(^M,Q)-?,k)\ M(^M,o) )[ M(M,i) / 

where the sum is over Mb, Md, M^m,o)j ^(m,!), k and the constraints 2. and 3. 
apply. We use the integral representation 

2m J ^ z^+^ ~ ^ J\A-2n 

1 f (l + z)^ (A 

az 



27TiJ z^+^ \B 
to reduce (|58|) to 



l dzidz2Si{zi, Z2)S2{zi, Z2) ^^ ^mIi (59) 
2m J Z2Z1 
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where 

Convergence of all series is guaranteed provided that we take \zi\ ~ e°, \z2\ ~ e 
with e <^ 1 and (3/2 < a < 3/5/4. Then one finds 

^n^i^^2) - -4 -3 ^-2^2 

Zi Z2 Z,ZiZ2 



/3 



5'2(^1,^2) 



Z2(l-^l)'(l + ^l 



^2 ^1^2 ~ 2^2^! — Zi — zf 

The integration over Z2 is performed first. Only the simple pole of 5*2 (-Zi, Z2) 

Z2 



1 - 2zi 

is encircled by the contour, while the poles of 5*1 lie outside the path of integration 
if e is sufficiently small. Finally, after some algebra, (p9|) reduces to 



1 f^ si+zr-'ii-3z) 



2m J z^+'^{l-4:z) 
The last equality has been obtained by deforming the contour of integration around 
the simple pole at z = 1/4. The proof given here is just the upgrading of an old 
method [|18]. Unfortunately we haven't been able to extend the proof to > 5. 
What seems to make higher N harder is not so much the constraint ( |50| ) , which can 
be taken care of by a double contour integration like for = 5, as the complicated 
form of vac{n) and vac{M,m). 

Clearly, the content rule (|50| ) has some consequences on the nature of the exci- 
tation spectrum in the thermodynamic limit, as computed in the previous section. 
In the FM case, it has been shown that only and ec are nonvanishing, but (|50D 
entails that roots other than (a) can be added only by increasing Ma and with it 
the energy. Hence, the ground state is not degenerate, and moreover no one particle 
excitation can exist (at least in this sector). 



As to the AFM case, we solve (M) and (EDI) for Ma, Mi 



NM NMc {N - 2)M, 



N-3 p-1 



Ma = — ^ - ^ -Y^nMn-Y. (2^^ + 1)M(M,™) (60) 



n—2 m=0 

even 



(N-2)M (iV-2)M, NMd ^ , 

^ ^ ^ n = 2 m=0 

We restrict ourselves to the simpler situation of even number of lattice sites. Then 
( |60| - |6T|) show that multiplets and ra-strings, n < N ~ 1, can appear as one-particle 



excitations, whereas Mc + M^ is bound to be even. The ground state has Ma = 



Mh = which exactly fill the available vacancies (|5T|). Since the distributions 



36|) are even in A, the only contribution to the ground state momentum is 

, ^,n^ ^NM 
Pa = 7rMi°) = 



2 

As to the excited states, AP picks a contribution (M^ — Mj^^^)7r. Owing to the 
pairing (|43| - ^5D and the arbitrariness of 27r in the definition of momentum, we have 



the following dispersion relations for the elementary excitation 
1. "Bound state" (c) plus hole in (a) 

^(A) = 2ec(A) P(A) = 27r,(A) + ^ P G (0, Nn) 



2. "Bound state" (d) plus hole in (a) 

E{X) = 2erf(A) P(A) = 27r,(A) + ill^^ p e (0, {N - 2)7r) 



3. (n, -|-)-string paired with (ra, — )-string, n = 2, 4, . . . — 3 

E{\) = 2e(„,±) (A) P(A) = 27r(„,±) (A) + mi P e (0, 2n7r) 

4. Multiplets m = 0, 1, . . . p - 1 

^(A) = e(M,™)(A) P(A) = 7r(Af,™)(A) + 7r(2m + 1) P G (0, (4m + 2)7r) 
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Calling P the upper limit of the momentum, we have in all four cases 

N 



E 



-P P^Q- 



and 



E 



N-1 

so that the group velocity is always \v 



N 



N — — 

(P-P) P^P 



N 
N-1 



5 Thermodynamics and central charge 

In the thermodynamic limit, the states of the spin chain are described by the density 
of rapidities pj and the density of holes p^^^ (alternatively the density of vacancies 



(Tj = Pj + p^ 



related by 



277 



Pk 



(62) 



The factor (±)''(j) has been included because Z{\) can be increasing or decreasing, 
so r(a) = r{b) = and r(j) = 1 if j 7^ a,b. It is well known |T0| , p!5| , p!9[| that, 
at finite temperature, the equilibrium state is determined by a system of nonlinear 
integral equations obtained by minimizing the free energy functional. The standard 
method, though, has to be suitably generalized to the present situation because not 
all densities pj, pf'' are independent. The following discussion is closely patterned 
after where the = 3 case has been dealt with. 

We assume that, not only for low lying states over the FM and AFM vacua, but 
for all states of the chain, the following constraints hold 



Pc = P^a^ and Pd = p^P 
pa = pf^ and Pb = p^d^ 

P{n,+) = Pin-) = Pn 



(63) 
(64) 
(65) 



Notice that (|63D and (|6J) imply cTc = aa and ad = crb- It is convenient to treat FM 
and AFM separately P]. We begin with FM and eliminate pc, Pd from the equations 
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with j 7^ c, d. We define the Fourier transform by 



and introduce the compact notation [n] = sinh(^) and {n} = cosh(^). We then 
find 

= "'j ~Y1 "^J,k * Pk (66) 

k 

where J2' is defined to extend to all types of roots other than c, d and 

= ~ 2{l}[iV- 1] = ^' ^ k = a,b 

[n] 

Tj,n{q) = -Tn,j{q) = _ j = a,b 

Tj^(M,m){q> = -T(M,m),j[q> = _ J = «, 

^ 2{l}[min(n,nO][A^- 1 - max(n,nO] 

T„x(?) = [W^] ^ ^ 

f_(g) = [N-n~ l][n -l] + [n][N-n- 2] 



,m),n 



[1][N-1] 

2{l}[min(ra,2m + l)][iV - 1 - max(n, 2m + 1)] 



im-i] 

2{l}[2min(m,m') + 1][A^ - 2 - 2max(m,m')] 
[2m][N - 2m - 2] + [2m + 1][N - 2m - 3] 



J-{M,m),{M,m')\q) = TT^t T\ [m TJl ) 

,m) 



[1][N-1] 

Note that the parities of the n-strings have disappeared because we find <J(n,+) = 
(T(n_) that implies, with (^), P(^''+) = p[n~)^ ^^^^ °^^y °f parities 
is independent and needs to be kept. Instead, replacing (|^) in (|62D with j = c,d 
yields Uc = cTq and ad = (Tb, showing that (0) actually follows from ( |62D and (|63|) . 
The energy density functional 
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is cast in the effective form by ehminating pc and pd by means of ( |5BD and then using 

E{p)=7rY.'J dXpj{X)aj{X) + eo (67) 

j 

where, hke in (p^), J^'j extends to (a), (6), (n), (M, m) and cq is the ground state 
energy density. Note that the functions aj (actually, in the case at hand, only one 
is nonvanishing) are closely related to the dressed energies and that, even in the 
effective form (^), the energy functional depends on pj only and not on p^'^'*. This 
is the reason why we have chosen to eliminate Pc and p^. In the entropy density 
functional |Ty], the degrees of freedom are also reduced by ( p^ - p5| ), so the effective 
expression is 

S{p) = J2' J K' log - Pj log Pj - pf^ log pf^] (68) 
j 

The minimum condition for the free energy density functional F{p) = E{p) — TS{p) 
is the system of integral equations 

ej = Tiaj + 'Tkj * log(l + e"^*^/^) (69) 

k 

where e~^/"^ = p/p^^\ The free energy itself is 

F = -TY.' j d\aj{\)\og{l + e-'^'^){\)+eo (70) 

It can be seen that the finite temperature dressed energies ej as defined by (|69D 
reduce, at T = to the actual energies of the excitations, rather than to the dressed 



energies as defined by (|29D . This is because the pairing "a-hole in c" has been 
included in the calculation imposing the constraint (|63| - |6^ ). We are particularly 
interested in the behavior of the specific heat at T 0. The regions at |A| ~ 
— 2^1ogT give the leading contribution to the T ^ 1 asymptotics. We define 
shifted functions /*(A) = /(A — ^^logT) and compare (|66|) and (|69D in the limit 
T ^ 0, finding 

pr-(-r»«^(i-/(|))i m 
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where f{x) = (1 + e^) ^. These relations can be inserted in the entropy, so that the 
leading term at T — can be written solely in terms of scaled dressed energies 0(A) 

S = Y: '^i-P'^ /7 ^^^t^^*^^) + + - + e-^O] (73) 

0(A) = limle(A-i^logT) (74) 



From their definition and (|69D, the scaled dressed energies solve 

0,(A) = 4iVe-2^^5,- , + %j * log(l + e-'^^)(A) (75) 

k 

We perform a change of integration variables in ([73|). When j = a,b —>■ g{(l)j) = 
(1 + e'^^y\ while for j ^ «, & 0i ^ /(0i)- Then 
2T 

^ = -]^{ E [^(^7(0j(+oo)))-L(^(0,(-oo)))]+ 5: '[L(/(0,(+oo)))-L(/(0,(-oo)))]} 
where i^(x) is the dilogarithmic Rogers function 

L(x)=-i rrf,(^+^°^(^-^)^ 



2 Jo 1 - y y 
The limiting values 0(±oo) can be obtained from the system of nonlinear (ordinary) 
equations to which ( |75|) reduces when A ±oo. The solutions are, at A — >^ +oo 



M+oo) = M+oo) = - log{N - 1) 
0„(+oo) = log((n+ 1)2-1) 
0(M,m)(+oo) = log((2m + 2)2 - 1) 



and at A 



-oo 



-oo 



-oo 



(M,m) \ 



+ 00 



log( 



sm 



• 2/(n+lV^ 

log(^^_\^ ' 



sm 



-oo 



log( 



sm 



2/ (2m+2)7r- 



sm 



N+2 



sm 



' 2n ' 
^ Af+2' 



^ Af+2' 



-11 



1) 



23 



By means of the identities 

f i((^) + 2L(i) = L(l) = ^ (76) 
fLf "°''^' l 2(JV-1) 

the leading term of the specific heat C = turns out to be 

from which we find the central charge (|^. 

The AFM case differs in the fact that, in order to avoid having hole densities 
in the effective energy functional, we prefer to eliminate pa and from (|62|). The 
result has the same form as (^) but now ^ is symmetric and 



[2]{N-l} "^^^ [2]{iV-l} 
'^nWJ — rn r AT TT '^{M,m) 



[l]{iV-l} 

[N-2] 



Tc,M) = TcAq) = Td,d{q) 



2[l]{iV-l} 



-Lc,n{q) = J^d,nW = JTYTTT — c,{M,m){(lj = d,{M,m){(l) 



Tn,n{q) 
Tn,{M;m) 



[l]{iV-l} ^'^'^-'-^^-^ -a,,..,m,,-,j [l]{iV-l} 

n')]{l}{A^- 1 -max(ra,n')} 
n]{N -n-2} + [n- 1\{N - n - 1} 



rf 2[min(n,nO]{l}{A^- 1 - max(n,nO} 
7^nx(^) = [l]{iV-l} ^^^^ 



[l]{iV-l} 

2[min(ra, 2m + l)]{l}{iV - 1 - max(n, 2m + 1)} 



^m) ,{M ,m) 



[1]{N-1} 

2[2min(m,m') + 1]{1}{A^ - 2 - 2max(m,m')} 
[1]{N-1} 

[2m]{N - 2m - 2} + [2m + 1]{N - 2m - 3} 



J-{M,m),{M,m')W = FTTTTT T^ [m ^ m ) 



[l]{iV-l} 



As to the equations with aa and at in (|6^) , once (|6^) is replaced we find cTq = (Xc 
and (Tfo = o"d, i.e. either one of (|63| ) and (Q) implies the other when inserted in (|6^). 
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With the new range of J2' and the new definitions of aj and Tj^k, (p7|), (pSD, (|5^), 
and (^) are formally the same (of course cq in (^) is now the AFM ground state 
energy density, as given in ( p7| ) ) . Again T = in (^) yields the energies of the zero 
temperature excitations, i.e. the zero temperature dressed energies after they have 
been paired by the correlations 



In the T <^ 1 limit, the relevant region is |A| ~ — ^^^j^^logT, so the shifted 
functions will read /*(A) = /(A — logT) and ( [n] - |7^) is replaced by (/(x) is still 
(l + e")~i 

(78) 

Defining scaled dressed energies 0(A) as in (^), we see that they take on the values 

at A +00 

0e(+oo) = M+oo) = log{N - 1) 
0„(+oo) = log((n + 1)2-1) 
0(M,m)(+oo) = log((2m + 2)2 - 1) 

and at A ^ — oo 

0j( — OO) = +00 Vj 

Using ( [78| - |79D , the entropy is now reduced to 

s = T.' J^y^^m^) iog(/(0.)) + (1 - /(0.))iog(i - fim 



2T(N - 1) 

^ ^.E'[L(/(0,(+oo)))-L(/(0(-oo)))] 



Nir 

3 



and the identity (^) is sufficient to determine the specific heat leading term 

7rr(iV-l) 



c 



3N 

N 



from which, since the velocity of the excitations is ^^/zj' c = 1 follows. 
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Appendix A 



In this appendix we prove that the rules ([S^ - ISq) imply the correlations of rapidities 



at least for finite energy excited states over FM and AFM vacua. We begin with 
FM and consider states having an arbitrarily high but finite number of holes in the 
(c) distribution and an arbitrarily high but finite number of other roots. We then 
have, neglecting terms of order 0(1/M) 

Z,(A) = ;it,(A)-^e,,,*pW(A) (80) 
^a(A) = ^aA)-i^e,,,*pf(A) (81) 

It's easy to see, from their Fourier transform, that Z'^{X) = —Z'^{X) and therefore 

Z,(A) = -Za{X) 

since they are both odd in A. Furthermore, ( PU] ) and (|HT| ) imply that they are both 
monotonic. But, by definition, 

J- a _ . . ih J- a 



SO from (^) we have 



Z,{\^;^) = -Z,(a|?'^)) = Z„(A(?'^)) (82) 



and since Za is monotonic, (|32|) follows. 
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In the AFM case, one considers, again discarding terms 0(1/M), 

= ^t,{x) - i^e,,, * pW(A) - i^e,,, * pr(A) 

with j = a, b, (n, +), (n, — ). In this approximation 

Za{X) = -Z,{X) 

Z,{X) = -Zd{X) 

^(n,+)(A) = ^(„_)(A) 

and all Z-functions are monotonic, so that, arguing as in (^), the rules (|5^ - |56D 



imply the pairings (^3|-^ 



Appendix B 



Let's suppose that the content rule (|50D and assumption (|48|) hold. We want to 



prove that the number of vacancies as defined by (|47| ) coincide with the more precise 



definition. Using (^q ) and (|50|) we see that the oddness table (p^ is modified 

j{a)j{c) ^ integer(half - odd) if M + + Ma = odd(even) 
j{b)^ j{d) ^ integer (half - odd) if Mb + Md = odd(even) 
j(n,+)^ = integer (half - odd) if M„ = odd(even) 

Furthermore, with (|46|) and (^O]) we find 

M + M, + M, M + Mb + M, 
MZ„(+oo) = ^ MZc(+oo) = ^ 

MZbi+oo) = ^ MZdi+oo) = ^^-^ 

Consider the case M + Mb + Md =odd. Then l'-"^ are integers and 

, _ + + ,83) 

It is easy to see that the number of integers in [—1^1^, Imlx] is M + Mb + Md- On 
the other hand, if M + M^ + Md = even, /'^"^ are half-odd. Again (|83|) holds, and 
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again, the number of half-odd numbers in [—1^1^, Imlx] is M + Mb + Ma- This is 

exactly the number of vacancies found from M{Za{+oo) — Za{—oo)). The same 

argument applies, with some obvious changes, to Zt,, Z^, Z^, noticing that the last 

two are decreasing and Imax is found from Z{—oo) = —Z{+oo). 

As to n-strings, the result is the same for {n, +) and {n, — ) so we consider only 

one. We find that the limit at A ^ — oo can be expressed 

M7 ( \ n(M - Me - Mg) ■ ( > (2n - 1)M„ 
MZ„(-oo) = 2^ mm(n ,n)M„/ 

- min(n, 2m + l)M(M,m) = 4"'' + 

where /q"'' is an integer. Again, we inspect the two possible cases. If M„ is even, 
1^^^ is half-odd and 

41 = + ^ (84) 
and the number of vacancies in [—Imax, Imax] is 2/o"''+M„, the same as M(Z„(— oo) — 
Zn{+oo)). Likewise, when Mn is odd, 1*^"^ are integers and (|^) still holds, giving 
the same number of vacancies. 

Finally, with multiplets, the Z-function is decreasing and from its definition, (^61) 
and (|50|) we have 

MZ(M,m)(-oo) = (2m + l)(M-Mc-Mrf)-2^min(n,2m+l)M„ 

n 

- J2 (4min(m, m') + 2)M^M,m') - (4m + l)M^M,m) (85) 

From (pop and (^Tj) we see that M — — is necessarily an even integer, so we 
can rewrite (^) as 

MZi^M,m){-00) = 24*'^'™^ - M(M,m) 

j{M,m) i^gjj^g integer. As shown in (p2D , J^^-'^'™-) must always be integers, still in the 
case = 3 ||13|, where only pairs (m = 0) are present, it was observed that 
is even(odd) when M(^Mfi) is odd(even). Let's assume this restriction hold for any 
A^ and any m, and inspect the two possibilities. If M(^M,m) is even, clearly 

= 2/o - M(M,™) - 1 (86) 
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Counting the half-odd numbers in [-4^^^ 4^^^] we find 24^^'""^ - M(M,m). If 
M(M,m) is odd, (|86|) still holds, but now /^i™-* is even and the number of even 
integers in the interval is the same. We conclude 

M 

Vac{M,m) = 2/o - M^M,m) = Y(^{M,m)(-00) - Z(M,m)( + 00)) 



which accounts for the factor 1/2 used in (51) 
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